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ABSTRACT 

Observations of the Cosmic Microwave Background (CMB) have revealed an unexpected quadrupole-octopole alignment along a preferred 
axis pointing toward the Virgo cluster. We here investigate whether this feature can be explained in the framework of the concordance model by 
secondary anisotropies produced by the non-linear evolution of the gravitational potential, the so-called Rees-Sciama (RS) effect. We focus on 
the effect caused by the local superclusters, which we calculate using a constrained high-resolution hydrodynamical simulation, based on the 
IRAS 1.2-Jy all-sky galaxy redshift survey, which reproduces the main structures of our Universe out to a distance of 1 10 Mpc from our Galaxy. 
The resulting RS effect peaks at low multipoles and has a minimum/maximum amplitude of -6.6/jK/1.9/jK. Even though its quadrupole is 
well aligned with the one measured for the CMB, its amplitude is not sufficient to explain the observed magnitude of the quadrupole/octopole 
alignment. In addition, we analyze the WMAP-3 data with a linear matched filter in an attempt to determine an upper limit for the RS signal 
amplitude on large scales. We found that it is possible to infer a weak upper limit of 30yuK for its maximum amplitude. 
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1. Introduction 

The Wilkinson Microwave Anisotropy Probe (WMAP) ob- 
servations of the Cosmic Microwave B ackground (CMB) 
dSpergel et all l2003t. Iffinshaw etail l200fjh turned out to be 
in very good agreement with the predictions of the concor- 
dance A-CDM cosmology. However, the WMAP data also 
present a number of unexpected features at large angular 
scales. In pa rticular, these i nclud e the small amplitude of the 



quadrupole dSpergel et aU I2003I) . as already no t iced t o be 



present within the COBE data by Hinshaw et al. ( 1996b . the 



alignment between the CMB quadrupole and o ctopole, which 



are se en to point toward the direction of Virg o (Tegmark et al 
2003), the planarity of the octopole dde Oliveira-Costa et al 
20041) . the alignment between the combined quadrupole and 
octopole mom e nta with the equinox a nd the ecliptic plane 



( Schwarz et al. . 2004 : Copi et al. . 2006h . the localized source 
of non-gaussianity (the cold spot) o n angular scale ~ 10° 



( Viel va et all l2004t ICruz et all |2005), and the asymmetry of 



the large-scale power between the two galactic hemispheres 



(Eri ksen et all 2004; H ansen et all 120041) . The latter seems to 



be the most relevant one. 

Different possible explanations for some of these features 
have been proposed, including a small u niverse with a cut-of f 
scale below the cosmological horizon dSpergel et al .1 120031) 



anisotropic universes such as Bianchi m odels dBunn et al 



1995 Kogut et al., ll997hlJaffe et all 2005), quin tessence mod 
els ( Dedeo et all 120031) or local dust-filled voids dlnoue & Silk. 
20061) . 



In the framework of the concordance A-CDM cosmology, 
we explore the possibility of explaining the above "anomalies" 
with the secondary CMB anisotropies produced by the evolu- 
tion of the gravitational potentials of the local cosmological 
structures. These anisotropies are usua lly split into the inte- 
grated Sachs- Wolfe (IS W) effect (ISW; ISachs & Wolfl 1 19671: 
Hu & Sugiyama, 1994) and into the Rees- Sciama effect (RS; 
Rees & S ciamai 1 19681; iMartinez-Gonzalez et all Il990l: ISeliaki 
1 19961) . which are produced by the linear and the non-linear evo- 
lution of the potential, respectively. We focus our attention on 
the RS effect although our simulation also includes the local 
ISW effect. However, the former is dominant over the latter in 
the volume we consider. We ignore distant structures beyond 
1 10 Mpc, because their contributi on to the RS effect is negligi - 
ble at low multipoles (see for e.g. Seli al 1 1 9961: Icooravl l2002h . 



A first study of t he RS imprint f r om the local universe has 
been carried out by ICoorav &Setol |2005). They used a toy- 
model with two spherical mass concentrations of 8xlO 15 /!~'M 
and 6.6 x 10 15 /z _1 M o to represent the Great Attractor and 
the Shapley concentration, and estimated a RS effect with 
a max imum amplitude of ~ 0.5^K. Similarly, iRakic et al 



(2006) modeled the Great Attractor with a spherical over- 
density (Lamaitre-Tolman-Bondi model) with the Local Group 
falling into it. They estimate a RS signal which may amount to 
AT/T ~ 10~ 5 and conclude that the RS effect is not compatible 
with the features found in the WMAP data. 

We use a high-resolution hydrodynamical simulation 
whose initial conditions were constrained to reproduce the den- 
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Fig. 1. Full sky map of the Rees-Sciama effect produced by the local universe within a sphere of radius 1 10 Mpc centered on our 
Galaxy. This CMB secondary anisotropy is dominated by the gravitational collapse of the most massive structures present in our 
constrained high-resolution hydrodynamical simulation. These are the galaxy clusters of Virgo, A3637, Centaurus, Hydra and 
Perseus. The maximum signal strength corresponds to an induced temperature shift of -6.6 fiK in the CMB. 



sity and velocity fields inferred with the IRAS 1.2-Jy all-sky 
galaxy redshift survey. The simulation gives a realistic picture 
of the local universe within a sphere of radius 110 Mpc cen- 
tered on our Galaxy. We use this simulation to obtain the ex- 
pected RS effect on large scales in the full sky, and to define an 
optimal linear filter for measuring it in the WMAP-3 data. 

The structure of the paper is as follows. In Section 2, we 
present the numerical simulation and the computation of the 
full sky map of the RS effect. In Section 3 we describe our 
analysis of the RS signal, and in Section 4 we propose a method 
based on a matched optimal filter to determine an upper limit 
for the amplitude of the RS effect in real data. Finally, we 
present our conclusions in Section 5. 

2. Simulating the RS effect due to the local 
universe 

2.1. The numerical simulation 

The results presented in this paper have been obtained by using 
the z — output of a cosmological hydrodynamical simula- 
tion of the local universe. We used initial conditions similar to 



those adopted bv lMathis et al.l d2002l) in their study (based on a 
pure N-body simulation) of structure formation in the local uni- 
verse. The galaxy distribution in the IRAS 1 .2-Jy galaxy survey 
is first smoothed with a Gaussian filter on a scale of 7 Mpc and 
then linearly evolved back in time up to z — 50, following the 



method proposed by Kolatt et al. I dl996l). The resulting field is 
then used as a Gaussian constraint dHoffman & Ribald, Il99ll) 
for an otherwise random realization of a flat A-CDM model, for 
which we assume a present matter density parameter Oo = 0.3, 
a Hubble constant Hq = 70 km/s/Mpc and r.m.s. density fluc- 
tuations of cr 8 = 0.9 in spheres of comoving radius 8 /i~'Mpc, 
with h = 0.7. 



The volume that is constrained by the observational data 
covers a sphere of radius ~ 110 Mpc centered on the Milky 
Way. This region is sampled with more than 50 million high- 
resolution dark matter particles and is embedded in a periodic 
box of size ~ 343 Mpc on a side. The region outside the con- 
strained volume is filled with nearly 7 million low-resolution 
dark matter particles, allowing a good coverage of long-range 
gravitational tidal forces that affect the high resolution region. 



The analysis by M afhis et al.l d2002l) demonstrated that the 
evolved state of these initial conditions provides a good match 
to the density and velocity fields of the large-scale structure ob- 
ser ved in the lo c al uni verse. Unlike in the original simulation 
by iMathis et al.l d2002l) . we however also follow the evolution 
of the gas distribution. To accomplish this we extended the ini- 
tial conditions by splitting each of the original high-resolution 
particles into a gas and a dark matter particle, having masses 
of 0.48 x 10 9 M and 3.1 x 10 9 M , respectively; this corre- 
sponds to a cosmological baryon fraction of 13 per cent. The 
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total number of particles within the simulation is then slightly 
more than 108 million. 

Our simulation run has been carried out with GADGET-2 



(Springel, 2005), a ne w version of the parall el Tree-SPH sim- 
ulation code GADGET dSpringel et all l200lh with an entropy - 
conserving formulation of SPH (Springel & Hernquist, 2002). 
The simulation neglected radiative cooling processes and em- 
ployed a comoving gravitational force resolution (i.e. the co- 
moving softening length) of 14 kpc (Plummer-equivalent), 
which is comparable to the inter-particle separation reached by 
the SPH particles in the dense centers of our simulated galaxy 
clusters. Previously, the same simulation has been used to study 
the propagation of cosmic rays and t o predict the SZ-eff ect 
from diffuse hot gas in the local universe [Dolag et aL ( 2005 allb ). 

In this work we neglect the contribution of all structures 
placed at a distance larger than HOMpc, because the vol- 
ume covered by our simulation is sufficiently large to account 
for all relevant sources of the RS effect that co ntribute to the 
low mu l tipoles we are interested in (see for e.g. ISeljakt Il996t 
Cooravl [2002, where the total RS power spectrum is shown 
to peak at 30 < I ). For example, the RS imprint of a galaxy 
cluster placed just outside the considered volume would sub- 
tend RS features with scales typically smaller than ~ 0.5°, i.e. 
corresponding to m ultipoles larger than I ~ 650 (see for e.g. 



Maturi et al 



2007). We note that the RS effect of galaxy clus- 



ters c an also be related to their gravitational lensing deflection 
field (Bi rkinshaw & Gulllll983l) whose typical scale is at most 
of the order of 1 Mpc. The inclusion of the neglected large vol- 
ume beyond 1 10 Mpc would be possible through la rger numer- 



ical simulations (see for e.g. ISpringel et all 120051) or through 
analytical modeling. However, a direct simulation approach 
would in general fail to reproduce the real matter distribution 
and dynamics, and hence fail to match with its induced RS ef- 
fect the actual orientation of the low multipoles of the CMB. 

2.2. RS full sky map 

Photons change their energy when they traverse evolv- 
ing gravitational potentials. This effect introduces secondary 
anisotropies in the CMB, which are usua lly split into two terms : 
the integrated Sachs-Wolfe effect QSWlSachs & WolfeL Il967l) 



and the Rees-Sciama effect (RS, Ree s & Sciamal 1 19681) . pro- 
duced by the linear and the non-linear evolution of the poten- 
tial, respectively. The resulting anisotropies can be calculated 

as 



AT 



<0) 



I 



<t>(0r, r/c) dr . 



(D 



where <fc refers to the time derivative of the gravitational po- 
tential, 6 is the position on the sky, r is the space coordinate 
along the line-of-sight, and r/c is the look-back time. These 
secondary temperature anisotropies provide precious informa- 
tion about the dynamical state and the evolution of the cosmic 
structures between the last scattering surface and us. In par- 
ticular, the signal produced by the local universe could con- 
tribute to the low multipoles of the observed CMB anisotropies 
in such a way to partially explain the deviations from the ex- 
pected Gaussian fluctuations as seen in the WMAP data. 



In order to estimate the secondary anisotropies due to the 
local universe from our simulation we consider a Newtonian 
approximation where the gravitational potential is defined as 



<D>(r, f) 



\r-r'\ 



so that its time derivative is 



<D(r, f) 



' J d\'p(r\ 



t)v(r',t) 



r — r 

r - r'\ 



(2) 



(3) 



where p and v are the matter density and velocity, respectively. 

We directly solve equatio n (O with a m odified version of 
the parallel code GADGET-2 dSpringell 120051) . whose original 
version solves numerically a similar integral to compute the 
gravitational forces. The RS map is finally obtained by adopt- 



ing the HEALPix pixelization of the sky dGorski et all 12005) 



with parameter n s id e = 512, i.e. the pixels on the sky extend 
over ~ 6.8' on a side, and by evaluating equation (Q]) at each 
pixel position. We compute the integral along the line-of-sight 
as a direct sum by sampling <j at regular intervals of 140 /z'kpc 
between the radius r m ; n = 1 Mpc and r max = 1 10 Mpc. In doing 
so we include only the high resolution region of the numerical 
simulation in order to avoid any border effects. In this volume, 
the resolution is good enough to properly sample the RS signal 
down to the typical cluster core scale. 

The resulting full-sky RS map due to the local universe is 
shown in Figure [TJ The main features are minima centered on 
those structures of the numerical simulation which correspond 
to the galaxy clusters of A3637, Virgo, Centaurus, Hydra and 
Perseus. The strongest of these signatures is related to Perseus 
(on the galactic equator, left hand side) with a negative peak 
of -6.6 //K, while the structure corresponding to the cluster of 
Coma shows a very small contribution to the overall RS signal. 
The RS effects of other prominent clusters are visible, but their 
overall impact is negligible. The RS map slightly favors nega- 
tive values because we sample a region of the universe popu- 
lated by forming structures with an ongoing gravitational col- 
lapse. 

3. Power spectrum analysis of the RS map 

The power spectra of our RS map, of the primary CMB, and of 
the instrumental noise of WMAP and of the upcoming Planck 
experiment are shown in Figure [2] The RS signal is dominated 
by large-scale features which are only above the instrumental 
noise level for a few multipoles (/ < 6 for WMAP and / < 30 
for Planck). However, the RS signal is always well below the 
CMB primary fluctuations even on the largest scales where our 
covered volume should be sufficient to give the full expected 
signal. A detailed comparison between the amplitudes of the 
dipole, quadrupole and octopole of the RS map and of the 



CMB , as observed in the 3-years WMAP data dHinshaw et al 



2006), is given in Table [TJ It is curious to notice the similar 
alignment of the quadrupole and octopole components of the 
RS and of the CMB, see Figure[3] But in any case, the RS effect 
from the local universe is far from being a relevant contribution 
to the CMB anisotropies because of its small amplitude, which 
corresponds only to a contribution at the percent level. 
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Fig. 2. Angular power spectra of the CMB, the RS effect 
from the local universe and the noise of WMAP (20 /jK and 
FWHM = 30') and of Planck (10yuK and FWHM = 5'). The 
RS effect is well below the amplitude of the primary CMB, but 
for a few low multipoles it is larger than the instrumental noise 
level. 

Table 1. Amplitude of the multipoles I = 1, 1 = 2 and I = 3 of 
the RS effect as derived in our simul ation, and of the CM B as 

E006h . 



measured by the 3 year WMAP data dHinshaw et al. 





WMAP 


RS 


1=1 dipole 


removed 


Ar, =] = 0.72 yuK 


1=2 quadrupole 


AT /=2 = 14.5^ 5 ytiK 


AT, =2 = 0.74 yuK 


1=3 octopole 


AT /=3 = 32+f yuK 


A7z =3 =0.51yt/K 



We conclude that the RS effect due to the local universe 
does not explain the deviations of the CMB from the expected 
primary anisotropies predicted by the concordance A-CDM 
model, i.e. the non-Gaussian features and the quadrupole- 
octopole alignment. However, we also found that the RS ampli- 
tude lies above the instrumental noise level for low I, and thus it 
may be possible to determine an upper limit for it, even if very 
weak, thanks to the present or upcoming full sky CMB exper- 
iments. On one hand a detection would help to understand the 
local cosmological structures and the details of their formation. 
On the other hand, if it was impossible to achieve a detection 
with an optimally designed filter then this would effectively en- 
sure that the RS effect from the local universe constitutes a neg- 
ligible foreground for any kind of data analysis of the CMB 
experiments. In the next Section, we therefore explore the fea- 
sibility of measuring the RS signal from CMB data by using a 
linear matched filter defined on the sphere. 

4. Measuring the RS effect 

We wish to construct a linear matched filter allowing the RS 
effect from the local universe to be extracted from millimetric 
observations of the CMB. We want the filter to suppress, as 
well as possible, the primary anisotropies of the CMB as well 
as the instrumental noise. In this work we ignore other possible 
sources of noise such us the residuals of the galactic foreground 
separation. Because of the significant level of such contamina- 



tions of the CMB primary anisotropies we only expect to obtain 
an upper limit constraint. 

4.1. Non symmetric matched filter on the sphere 

The filter used in this work is based on a template for the 
expected signal and on the power spectr u m of the noise. 
It was proposed by Haehnelt & Tegmark ( 19961) in a flat 
sky approximation to estimate the peculiar velocities of 
galaxy clusters through their kinetic Sunyaev-Zel'dovich ef- 
fect. Subsequently, it has also been proposed as a tool to es- 
timate the mass of galaxy clusters th rough weak lensing of 
the CMB dSeljak & Zaldarriagai l2000|) and for cluster d etec- 
tion in optical weak lensing surveys dMaturi et al. . 2005). An 
extension of the filte r to the sphere was first introduced by 
Schaef er et al. (2004), assuming axially symmetric models for 
the expected signal. Since we are interested only in the ampli- 
tude of the RS signal from the local universe and not in recov- 
ering a map of the full sky, it is straightforward to generalize 
the filter for non-axially symmetric models, like for the RS ef- 
fect. Our approach hence avoids the significant difficulties in 
convolving the sky with a highly asymmetric filter. 

We model the measured data d(6), which contains a signal 
s(0) = A t(0) contaminated by some noise n(6), by 



d(6) = s{0) + n(0) = At(6) + n(6) 



(4) 



In our case, n(6) is mainly given by the sum of the CMB pri- 
mary anisotropies and the instrumental noise. Both can be suf- 
ficiently well approximated with Gaussian random fields with 
zero mean. The noise correlation function is 



(5) 



where the asterisk denotes the complex conjugate, 6 is the delta 
function, and C; = C/ cm b + C/ ; j nst is the total angular power 
spectrum of the noise. 

For convenience, all quantities are derived in the harmonic 
space. We wish to construct a linear filter V F(0) which yields an 
estimate A for the amplitude A of the signal, viz. 



A = 



+<x> I 

Z 

1=0 m=-l 



Z dlnX lm , 



(6) 



where d\ m and ^F/,,, are the harmonic coefficients of the data and 
of the filter function, respectively. The sum over I is usually 
restricted to the maximum multipole given by the resolution of 
the observations. 

The estimate is required to be unbiased and optimal, i.e. the 
quantity 



£ = (A - A) = A 



Z Z r "»^ 

1=0 m=-l 



(7) 



has to vanish and the r.m.s of the measurement 
o- 2 = ((A - A) 2 ) = ^ 2 + Z Z C/ mm|2 ' 

/=0 m=-l 



(8) 
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Fig. 3. Comparison between the quadrupole and the octopole of the CMB as measured by WMAP (top panels) and of the RS 
effect from the local universe (bottom panels). 



has to be minimal. The filter *F which satisfies these two condi- 
tions is the one which minimizes the action L — a 2 + Ab, where 
A is a Lagrange multiplier. We thus obtain 



"Vim = 



EE 

1=0 m=-l 



Ci 



Tim 
Q 



(9) 



This shows that the filter is most sensitive for those spatial fre- 
quencies where the signal r is large and the noise C/ is small. 
Note that this filter uses the full information available on the 
expected signal r and not only its power spectrum, as for ex- 
ample a Wiener filter would do. 

4.2. An upper limit for the RS amplitude in the 
WMAP-3 data 

Before applying our procedure to real data, we have tested the 
described filter on mock data where we combined the full sky 
RS map (see Section l2~2l from our simulation of the local uni- 
verse with a realization of the primar y CMB fluctuations com- 
puted with CMBEASY JPoranl l2003h for a standard WMAP-3 
cosmology. The resulting map was then degraded by adding an 
instrumental noise with the power spectrum 



C, 



/,inst 



w exp 



/(/+ 1)^WHM 



81n2 



(10) 



where w 1 := (<r r #fwhm) 2 dKnox , 1995 ). The resolution and 
sensitivity were fixed according to the characteristics of the 
WMAP satellite (6»fwhm = 30', = 20//K) and to the up- 
coming ESA space mission Planck (#fwhm = 5',cr r = 10//K). 
A convolution with a Gaussian kernel with the same FWHM 
is applied to the simulated data. The dipole component of the 
signal is finally removed because of the strong contamination 
which would be caused by the proper motion of the solar sys- 
tem. We ignore other contaminants such as residuals of the 
foreground subtraction or extragalactic point sources. 

The resulting map is finally analyzed with the filter de- 
scribed in Section l4~Tl For the filter template r we use the same 
RS map which was also included in the simulated observations. 
This is a strong simplification, but here we only aim at testing 
the performance of the filter when all the conditions for its 
application are optimally satisfied. The results are as follows: 



Expected signal 
WMAP (simulation) 
Planck (simulation) 



A = 1.9 juK 

A = 3.4 ± 13/iK 

A -2.2 + 12yuK 



It is thus only possible to infer a weak upper limit for the am- 
plitude of the RS effect due to the local universe, because of the 
strong contamination given by the CMB primary fluctuations. 

Finally, we apply the filter, as previously d escribed, to the 
3 -years ILC map of the WMAP experiment dHinshaw et al 
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Fig. 4. Filter function used to estimate the upper limit for the 
maximum amplitude of the RS effect produced by the local uni- 
verse. 



2006). In this case, since the ILC map is convolved with a 
Gaussian kernel with FWHM = 1°, we convolved the template 
t in the same way before the filter derivation. This convolution 
of the template has the implicit advantage of suppressing small 
features in the RS template which may not be present in the 
structures of the real universe. The resulting filter is shown in 
Figure |U 

To estimate the variance of the measurement, we refrain 
from using Equation ([8j because it supposes that the template 
is identical to the expected signal and thus Eqn. © would 
yield a lower limit for the noise amplitude. We rather estimate 
the noise level by applying the filter many times with the RS 
template randomly rotated on the unit sphere. While the mean 
of these randomized measurements should be consistent with 
zero because the rotated templates do not match the expected 
signal, the dispersion around zero provides a good estimate of 
the noise level. We obtain 



WMAP (ILC map) 



A = 3 ± 26 juK 



which implies an upper limit for the amplitude of the RS 
effect in the actual CMB data of A max = 30 /iK. This is not very 
restrictive compared to the expected signal strength of A exp = 
1.9/iK. 

5. Conclusions 

We estimated a realistic map of the Rees-Sciama effect causd 
by the evolution of the gravitational potential in the local uni- 
verse within a sphere with a radius of 1 10 Mpc centered on our 
Galaxy. To this end we used a high-resolution hydrodynami- 
cal simulation whose initial conditions were constrained to re- 
produce the density and velocity fields inferred from the IRAS 
1 .2-Jy all-sky galaxy redshift survey. 

We estimated a minimum/maximum amplitude of the lo- 
cal RS effect of -6.6 yuK/ 1 . 9 pK, which is a factor of 10 larger 
than the 0.5/iK predicted byE ooray & Se to (2005). The power 
spectrum of the RS signal lies above the instrumental noise of 
WMAP and Planck for low multipoles but is always well be- 
low the CMB primary anisotropies. Interestingly, the RS sig- 



nal's quadrupole and octopole moments are well aligned with 
those of the CMB primary fluctuations. However, the small 
amplitudes of these RS multipoles are not sufficient to ex- 
plain the quadrupole-octopole alignment observed by COBE 
and WMAP in the CMB. 

In addition, we attempted to measure the RS effect on large 
scales in mock CMB observations as well as in real data using 
a matched filter technique. We found that it is only possible to 
place a rather weak upper limit for its amplitude even under 
optimistic assumptions. In ana lyzing the 3-years ILC map of 
WMAP dHinshaw et all 120061) this upper limit turned out be 
Amax = 30 yuK, which is not very restrictive compared with the 
expected signal of A exp = 1.9 yuK. 
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